Lowest-lying Tetra-Quark Hadrons in Anisotropic Lattice QCD 
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We present a detailed study of lowest- lying q 2 q 2 hadrons in quenched improved anisotropic lattice 
QCD. Using the irn and diquark-antidiquark local and smeared operators, we attempt to isolate 
the signal for I{J P ) = 0(0+), 2(0+) and 1(1+) states in two flavour QCD. In the chiral limit of 
light-quark mass region, the lowest scalar Aq state is found to have a mass, m{^° = 927(12) MeV, 
which is slightly lower than the experimentally observed /o(980). The results from our variational 
analysis do not indicate a signature of a tetraquark resonance in I = 1 and 1 = 2 channels. After 
the chiral extrapolation the lowest 1(1 + ) state is found to have a mass, m^ 1 = 1358(28) MeV. 
We analysed the static 4q potential extracted form a tetraquark Wilson loop and illustrated the 
behaviour of the 4q state as a bound state, unbinding at some critical diquark separation. From our 
analysis we conclude that scalar 4q system appears as a two-pion scattering state and that there is 
no spatially- localised 4q state in the light-quark mass region. 

PACS numbers: ll.15.Ha, 12.38.Gc,11.15.Me 



I. INTRODUCTION 

The concept of multi-quark states has received re- 
vived interest due to the narrow resonances in the 
spectrum of states. The recent experimental observa- 
tion of several new particles as candidates of multi- 
quark hadrons are expected to reveal new aspects of 
hadron physics [l], 0, H, 3- Among these discover- 
ies, the tetraquark systems are also interesting in terms 
of their rich phenomenology, in particular the scalar 
mesons which still remain a most fascinating subject of 
research. The five ++ isoscalar mesons below 2GeV: 
/o(400 - 1200),/o(980),/o(1370),/ (1500) and / (1760) 
are subject of conflicting interpretations. /o(1370) is as- 
signed to be lowest qq meson, / (1500) and / (1760) 
are expected to be the lowest scalar glueball or an ss 
scalar meson. The differences in the interpretations man- 
ifest themselves in the flavour wave functions of scalar- 
isoscalar mesons. To resolve the interpretational dif- 
ferences would require a detailed analysis of J/yb de- 
cays into vector plus scalar mesons and the radiative 
decays like /o — ► and /o — > "f(f> that constrains the 
flavour structure. Such experiments are feasible at BESII 
[E Hi S Hi and BESIII with good chances for all the chan- 
nels in radiative J/ip decays. The tetraquark systems 
are also interesting in terms of the nearly stable scalar 
charm state Dg(1980), the narrow £>q(2450), the ana- 
logue of an(980) with open charm, AT(3872) and F(4260) 
[ill id [Til . Ila . HH . From the consideration of their masses, 
narrow decay width and decay pattern, these states can- 
not be regarded as the simple mesons of cc or cs but 
conjectured to be 4g or molecular states [I'll. [iH [iH] . 

While a phcnomcnological model may often be helpful 
in obtaining a qualitative understanding of the data, lat- 
tice QCD, at the present status of approximations, seems 
to provide a trustworthy guide into unknown territory 
in multi-quark hadron physics [l7l Il8j . It is expected 



to play an important role in revealing the real nature 
of multi-quark hadrons. There are several, mainly ex- 
ploratory, lattice studies of tetraquark systems available 
[H, HI ED, HI HI, HI [H, HI , with mixed results. Most 
of these calculations claim the existence of a bound scalar 
1 = tetraquark state whereas Ref. [26j] has observed no 
evidence for a \q resonance. One may naively interpret 
the negative results obtained in Ref. (2|| as a conse- 
quence of their combined analysis of maximum entropy 
method with hybrid boundary condition method on their 
choice of operators, which would strongly couple to two- 
pion scattering states. Thus, it is important to use var- 
ious interpolating operators and variational techniques 
on more accumulated data to make a precise prediction. 
Our approach will be to use a number of different in- 
terpolating fields and smeared operators to enhance the 
low-lying spectra. 

In this study, we provide a detailed analysis on 
tetraquark hadrons from the combination analysis of 
variational method and smearing on various interpo- 
lating fields on improved anisotropic lattices. Us- 
ing the quenched approximation, and discarding quark- 
antiquark annihilation diagrams, we construct q 2 q 2 
sources from multiple operators. We exclude the pro- 
cesses that mix qq and q 2 q 2 and allow the quark masses 
to vary from small to large values. In the absence of 
quark annihilation, we do not expect any mixing of q 2 q 2 
with pure glue. Thus we can express the q 2 q 2 correlation 
functions in terms of a basis determined by gluon and 
quark exchange diagrams only. In the process of search- 
ing for a bound state it is essential to explore a large num- 
ber of interpolating fields having the quantum numbers 
of the desired state. Explicitly, one needs to construct 
an interpolating field which has significant overlap with 
the 4g system. However, any (qqqq) operator must cou- 
ple to hadronic states with the same quantum numbers 
(for example, the 7T7t scattering state). It is necessary to 
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disentangle the lowest-lying 4q states from the tttt scat- 
tering states, as well as the excited Aq states. To this 
end, we adopt a variational method to compute a n x n 
correlation matrix from the interpolating fields and from 
its eigenvalues we extract the masses. 

The rest of the paper is organized as follows: in Sec. II, 
we discuss our choice of interpolating fields for the sup- 
posed scalar tctraquark state and outline our construc- 
tion of the correlation matrix analysis. The technical 
details of the lattice simulations as well as the actions 
used in this study are discussed in Sec. III. It is essential 
to use the improved action that displays nearly perfect 
scaling, since large scaling violations could lead to a false 
signature of attraction. We therefore perform our simula- 
tions on improved lattice actions with tadpole-improved 
renormalised parameters that provide continuum limit 
results at finite-lattice spacing. We present and discuss 
our numerical results in Sec. IV. Finally, we perform the 
study of the inter-quark potential in our 4q system to 
establish a connection between connected 4q state and 
the "two- meson" state around the cross-over. Sec. V is 
devoted to our summary and concluding remarks. 



II. LOWEST q 2 q 2 STATES ON THE LATTICE 

We propose interpolators designed to maximize the 
possibility to observe attraction between tetraquark con- 
stituents at relatively heavy quark masses. Two general 
types of operators are considered: those based on the tttt 
configuration, and those based on a diquark-antidiquark 
configuration. The simplest 7T7r-type interpolators for the 
1 = and 1 = 2 channels, respectively, are 

O{= (x) = [{(d a (x) l5 u a (x))(u b (x) l5 d b (x)) 

-(u <-> d,u <-> d)} + ^{(u a (x)~f 5 u a (x)) 

x(d b (x)j 5 d b (x)) ~(u^d,u^ d)}] , (1) 



and 



Oi= 2 {x) = (d a ix) 75 u a (x))(d b (x) 75 u b (x)). (2) 



The J p = 1+ operator with isospin I = 1 is evaluated 
from the pseudo-scalar and vector mesons fields and has 
the form: 



Ol=\x) 



\ [(d a (xHu a (x))(^(x)^ 5 u b (x)) 

-{d a (x) lb u a (x)){d % {x) li u h {x)) 
+ {d a {x) li u a {x)){d!'{x) li u b {x))\ . 



(3) 



The other type of interpolating field is one in which 
quarks and antiquarks are coupled into a set of diquark 
and antidiquark, respectively and has the form 



Oi(x) = e abc [u a CTd b }e dec [u d CTd, 



(4) 



Accounting for both colour and flavour antisymmetry, 
the possible Ts are restricted to within 75 and 7j. The 



operators O12 are similar to the naive pion (gi pion op- 
erator with the difference that in the former one has a 
contribution from gluon and quark exchange diagrams to 
the 7r — 7r four-point functions. The operator O4 is mo- 
tivated by the Jaffe-Wilzcek description of diquarks [27| , 
which is expected to have a small overlap with two-meson 
scattering states. 



A. Extraction of masses 

The mass of the ground state is extracted from the 
asymptotic behaviour of the two-point the temporal 4Q 
correlator 



C(t) = ^(O(f,t)Ot(0,0)), 



(5) 



where the total momentum of the 4g system is projected 
to be zero. To disentangle the lowest-lying Aq states 
from 7T7T scattering states this way, we adopt a variational 
method to compute a 2 x 2 correlation matrix from two 
different interpolating fields. 

We compute the propagators {{Oi) x {Oj) y ) with fixed 
y = (0, 0), and their time correlation functions, 



Ci J 'W = (X]tr[((O i )(x,t)O J )(0,0)) J 



(6) 



where the trace sums over the Dirac space, and the sub- 
scripts / and U denote fermionic average and gauge fi eld 
ensemble average, respectively. Following [l7ll28l.]29j we 
solve the eigenvalue equation 



C(t )v k (t ) = Xk(to)v k (t ) 



(7) 



to determine the eigenvectors Vfc(io) an d use these eigen- 
vectors to project the correlation matrices to the space 
corresponding to the n largest eigenvalues X n {to) 



h3 



1. 



(8) 



and solve the generalised eigenvalue equation for the pro- 
jected correlation matrix C™ . In practice, we extract our 
results from a 2 x 2 correlation matrix. The masses are 
extracted from the average of CV, by a single-exponential 
fit to Xi (t) for the range of t in which the effective mass 
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A(t) 



A(t + 1) 



(9) 



attains a plateau. However, to ensure the validity of our 
results, we compare them with those obtained using 



M eff = In 



\{t - 1) - \{t) 
A(t)-A(t + 1) 



(10) 
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B. Static 4q potential 



III. SIMULATION DETAILS 



The static tetraquark potential is extracted from the 
gauge invariant 4q Wilson loop. The SU(3) Wilson loop 
is constructed by creating a 4q state at a time t = 0, 
which is annihilated at a later time t. Following 0, [25[ 
we write the 4q Wilson loop as 



J_e a6c e de - f e a ' 6 ' c 'e d ' e '- f 'U'(x,x' ) l) aa '{7(x ) x',2 



/ n\bb' 



12 



x(7 G (x,y)^C/(y',y,3) d ' d C/(y',y,4) e ' e 
xU G ,(y',x'y' c ' (11) 
where the staple-like links C/(x, x',/c) are given by 



ig / dz^A^z) 



£/(x, x', k) — Pcxp 



where denotes the path from x to x' for quark line k. 
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FIG. 1: The Wilson loop for the tetraquark system. 

The junctions at x and y arc joined by a line-like con- 
tour 



tfc?(x,y) = Pexp 



ig I dz ■ A(z) 

G 



at t = 0. The tetraquark potential can be obtained 
from the large t behaviour of (Wig), where the ground- 
state contribution becomes dominant. To obtain the op- 
timal signal-to-noise ratio, we must suppress the contri- 
butions from the excited states, which may be done by 
APE smearing technique 30] . This involves replacing the 
space-like link variables Ui by smeared link variables. To 
tune the smearing parameters, we fix the smearing frac- 
tion to a — 0.7, as it is sufficient to fix a and tune the 
smearing level to uape- A typical value which enhances 
the ground-state component in the 4q Wilson loop was 
tiape = 8. The tetraquark potential is then extracted 
in the standard way from the large t behaviour of the 4g 
Wilson loop with a single-exponential fit: 

{W 4q ) = Ze- y ^ (12) 

in the range [t m i n , t max \. The effective potential is then 
extracted from 



V iq {t) = In 



W ig (t) 



W iq (t + 1) 



(13) 



In this way, we calculate the potential for our planar 
configurations. 



To examine q 2 q 2 in lattice QCD, we explore the im- 
proved actions on anisotropic lattices. With most of 
the finite-lattice spacing artifacts having been removed, 
one can use coarse lattices, with fewer sites and much 
less computational effort. It has been found that even 
on fairly coarse lattices these actions give good results 
for hadron masses by estimating the coefficients of the 
improvement terms using tadpole- improvement. Using 
a tadpole improved anisotropic gluon action [3l| . we 
generate quenched configurations on 16 3 x 64 lattice 
(with periodic boundary conditions in all directions) at 
(3 = 4.0. Gauge configurations are generated by a 5- 
hit pseudo heat bath update supplemented by four over- 
relaxation steps. These configurations are then fixed to 
the Coulomb gauge at every 500 sweeps. After discard- 
ing the initial sweeps, a total of 300 configurations are 
accumulated for measurements. Quark propagators are 
computed by using a tadpole- improved clover quark ac- 
tion on the anisotropic lattice [32|]. All the coefficients 
in the action are evaluated from a tree-level tadpole im- 
provement. We calculate the spectrum at quark masses 
as light as 0.004, which is much less than the hadronic 
scale, 1.0 GeV, associated with the scalar Aq system. 

On each gauge fixed configuration, we invert the quark 
matrix by the BiCGStab algorithm to obtain the quark 
propagator. Using anti-periodic boundary conditions in 
the time direction, we calculate the propagators on three 
source time slices on the same lattice. For the lattice 
analysed here, we place the wall sources on the first, sec- 
ond and third time slices. Since there is a possibility of 
correlation among the propagators with different source 
time slices, we average them and treat them as a single 
result in a jackknifc analysis. 

The long distance exponential falloff is determined by 
the lowest-energy hadron state to which the operator 
couples. For lattice QCD, one has options to use a lo- 
cal or nonlocal interpolating operator for hadrons, pro- 
vided its correlation function has a significant overlap 
with the state under consideration. To obtain a better 
overlap with the ground state, we used iterative smear- 
ing of gauge links and the application of the fuzzing tech- 
nique for the fcrmion fields [33| taking the physical size of 
the particle into account. The fuzzed quark field is con- 
structed to be symmetric in all space directions. Such a 
fuzzed quark field is only used at the sink, while the one 
at the source remains local. The application of fuzzing 
for two of the four quarks inside the q 2 q 2 flattens the cur- 
vature of the effective mass. The largest plateau in the 
region with small errors is obtained with fuzzed u- and d- 
quarks. We used this variant to calculate our correlation 
functions. 

We estimate the lattice spacing by linearly extrapolat- 
ing the p mass to the physical quark mass. The latter can 
be determined from the ratio of two non-strange hadrons. 
Since m u and md are fixed by the nucleon mass, in the 
non-strange sector, m p is a test both for finite-size effects 
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and quenching errors. The major source of discrepancy 
among the lattice spacings from different observables is 
the quenching effect. The obtained p and N masses are 
compared to the experimental values and show a devia- 
tion of less than 3 — 4% for the lattice size explored here. 
Such a variance can be considered as the usual quenching 
effect. Inspired by the good agreement of the p and N 
masses with the experimental values, the scale was set al- 
ternatively by the ratio ui^/uin. Using the experimental 
value 768 MeV for the rho mass, the spacing of our lattice 
is a s = 0.462(2) fm. The bare quark mass is determined 
by extracting the mass of the K meson. At m q at = 0.04, 
we obtained the mass of vector meson ttik = 553(2) MeV, 
which is in good agreement with the estimates obtained 
from other studies 0, [34|. Thus taking the strange 
quark bare mass to be m s at = 0.04, we choose the bare 
quark masses for u and d to be much smaller than m s , 
i.e., m q a t = 0.01,0.008,0.0065,0.0055,0.005,0.004. 
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FIG. 2: Effective mass of Aq in the 1 = (open triangles) and 
1 = 2 (open circles) channels at n t = 0.2635 (m q a t = 0.0065). 



IV. RESULTS AND DISCUSSION 

In this section we present the results of our lattice sim- 
ulation of tetraquark masses for the isospin I = 0, 1 and 
1 = 2 channels. In addition to extracting the masses of 
the q 2 q 2 states, we also study the mass differences be- 
tween the candidate tetraquark states and the free two- 
particle states. This analysis is actually important in de- 
termining the nature of the states observed on the lattice, 
and the identification of the true resonances. Hence a 
lattice signature which might be observed for tetraquark 
resonance is a negative mass splitting at quark masses 
near the physical regime. 

Fig. [2] shows the effective mass plot for the low- 
est scalar Aq (udud) for the typical hopping parameter 
n t = 0.2635. At sufficiently large t, contributions from 
excited states are diminished, and the correlator is dom- 
inated by a single state. Thus a plateau may serve as an 
indicator of the single-state saturation, hence allowing 
one to perform the single-exponential fit to the plateau 
region. In the region < t < 9, the effective mass for 
the 1 = channel decreases monotonically and reaches 
a stable plateau at 10 < t < 25. Beyond t ~ 25 the 
data become noisy. The fitting range [t m in , t m ax] for the 
final analysis is determined by fixing t max and finding 
the range of t m j„ where the ground-state mass is stable 
against t min . We choose one "best fit" which is insensi- 
tive to the fit range, has high confidence level and reason- 
able statistical errors. We then confirm this by looking 
at the plateau region of the correlator. Statistical er- 
rors of masses are estimated by the jackknife method. 
Since the analysed configurations were highly uncorre- 
cted, which we ensured by separating the analysed con- 
figurations by as many as 1000 sweeps, the statistical 
errors of our mass estimates are typically on the few per- 
cent level. With much of our data we find that the er- 
rors are purely statistical, and the goodness of the fit is 
gauged by the xV-^df- All fits are reasonable, having 



X 2 /Vdf < 1- For the 1 = 2 channel it was possible to 
find a fit region [£15,^35] in which a convincing plateau 
was observed. The effective mass is found to be stable 
using different values of t in Eq. @ , which suggests that 
the tetraquark ground state is correctly projected. 
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FIG. 3: Effective mass of I(J P ) 
0.2650 (m q a t = 0.005). 



1(1 + ) Aq state at K t 



For the I{J P ) = 1(1+) tetraquark state, Fig. El we 
observe that the exponential seems to oscillate in time at 
the lowest-quark mass at large t. One can still demon- 
strate the existence of a plateau in the effective mass 
and extract the mass with some reliability, although the 
errors are relatively large. Suppressing any data point 
which has an error larger than its mean value, a possible 
plateau is seen in the region 20 < t < 35 with reasonable 
errors, where the single-state dominance is expected to 
be achieved. To ensure the validity of our results, we 
compared them to those obtained using Eq. (|10[) . It was 
found that the evaluations of ([9]) and (fT0|) yield results 
very consistent within statistical errors. 

The chiral extrapolation to the physical limit is the 
next important issue. From the view point of chiral per- 



turbation theory, data points with smallest m\ should be 
used to capture the chiral log behaviour. Leinweber et al. 
[35l ] demonstrated that the chiral extrapolation method 
based upon finite-range regulator leads to extremely ac- 
curate values for the mass of the physical nucleon with 
systematic errors of less than one percent. We use a set 
of data points with smallest m\ to capture the chiral 
log behaviour. Fig. Q] collects and displays the resulting 
particle masses extrapolated to the physical quark mass 
value using linear and quadratic fits in m^, 



m h = a + bml, m h 



cm, 



The difference between these two extrapolations gives 
some information on the systematic uncertainties in the 
extrapolated quantities. Although, our quark masses are 
quite small, both linear and quadratic fits essentially gave 
the identical results. The chiral uncertainties in the phys- 
ical limit are significantly small (less than a percent) at 
our present statistics. 




FIG. 4: Chiral extrapolation of effective masses in the 1 = 
(open circles) and 1 = 2 (open triangles) channels. The solid 
curve denotes the two-pion threshold, and dashed curves are 
the linear fits to the data. 

At smaller quark masses, the data for the I — channel 
appear slightly below the two-pion threshold, by ~ 30 
McV, and seems to behave linearly in m\ rather than the 
two-pion threshold which shows a quadratic behaviour in 
m\ . The data for I = 2 lie above the two-pion threshold 
by ~ 40 — 50 MeV for all quark masses analysed here. 
If the slight difference from two-pion threshold can be 
explained by the the two-pion interaction, then the lowest 
scalar 4g state can be regarded as a two-pion scattering 
state. 

The results for the extracted mass for J p = 1 + are 
displayed in Fig. as a function of m\. The ground- 
state masses for the tetraquark and two-particle states 
are again very different. The mass of the tetraquark 
state extracted is consistently higher than the lowest two- 
particle state, namely 7r + p. This trend continues in the 
physical limit where the masses exhibit the opposite be- 
haviour to that which would be expected in the presence 
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FIG. 5: Chiral extrapolation of effective masses in I(J P ) = 
1(1 + ). The open triangles denote the Aq data and the solid 
diamonds denote the two-particle state (n + p). The dashed 
curve is a linear fit to the data. 



of binding. In the chiral limit the / = 
the two-pion threshold by ~ 70 MeV. 



1 state lies above 



TABLE I: The masses of the I(J P ) = 0(0+), 2(0+) and 
I(J P ) = 1(1 + ) 4g states in the lattice units for various values 
of Kt- 







J p 


= 0+ 


J p = l+ 




M p 


MZ=° 




Mlf 


0.2610 


0.850(2) 


0.892(4) 


0.956(6) 


1.305(4) 


0.2620 


0.741(2) 


0.815(5) 


0.877(6) 


1.174(6) 


0.2635 


0.687(3) 


0.758(5) 


0.808(7) 


1.094(7) 


0.2640 


0.606(3) 


0.709(6) 


0.756(7) 


1.012(7) 


0.2650 


0.572(3) 


0.663(7) 


0.712(6) 


0.949(8) 


0.2655 


0.529(4) 


0.627(8) 


0.668(8) 


0.884(9) 



Since quenched spectroscopy is quite reliable for the 
mass ratio of stable particles, it is physically even more 
motivated to extrapolate the mass ratio instead of the 
mass. This allows for the cancellation of systematic er- 
rors since the hadron states are generated from the same 
gauge field configurations and hence systematic errors 
are strongly correlated. The mass ratios at hand show 
a remarkably small scaling violation; -hence, we adopt 
an m^-linear extrapolation for the continuum limit. We 
also perform an m^-quadratic extrapolation to estimate 
systematic errors. Performing such extrapolations for 
all sets of masses, we adopt the choice which shows the 
smoothest scaling bahaviour for the final value, and we 
use others to estimate the systematic errors. 

The results of the mass difference analysis are shown 
in Table [IT] and illustrated in Fig. [6l The data behave 
almost linearly in m\ and both linear and quadratic fits 
essentially gave identical results. Again the contributions 
from the uncertainties due to chiral logarithms in the 
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FIG. 6: Extrapolation of the mass ratios of the 1 = (open 
circles) and 1 = 2 (open triangles) channels to the physical 
limit. The dashed lines are the linear fits in mt to the data. 



physical limit are seen to be significantly less dominant. 

The question whether the lowest-lying scalar Aq state 
extracted is a scattering state or a resonance is better 
resolved by analysing the ratio of the mass differences 
AM = M± q — 2M q q (between the candidate lowest scalar 
AQ state and the two-pion state) and the p mass. 

The non-zero m 2 . values of the ratio are within 0.01- 
0.02 standard deviations of the extrapolated zero pion 
mass result. This implies that the chiral uncertainties in 
the physical limit are less than 2%. The fitted results 
show that the mass difference AM and hence the mass 
ratio for the 1 = channel arc clearly negative and in- 
crease in magnitude as we approach the physical regime. 
This implies that the binding becomes stronger at light- 
quark masses with a general trend of negative binding as 
the zero quark mass limit is approached. In the phys- 
ical limit, the negative mass difference is ~ 25 — 30(2) 
MeV. Naively one may be tempted to interpret this neg- 
ative mass splitting near the physical regime as a lat- 
tice resonance signature for scalar Aq resonance with 
1 = 0. Using the chirally extrapolated value of the ratio 
m^/mp = 0.60(1), the 1 = state is found to have a 



mass of : 



-.1=0 
L 4q 



927(12) MeV in the physical limit. The 



lowest scalar Aq state appears to be slightly lighter than 
the experimentally observed /o(980). It still remains to 
verify whether analysis at relatively large quark masses 
would affect the manifestation of the bound state and aid 
to confirm the indication of the resonance. This analy- 
sis is discussed in the next section. On the other hand 
the 1 = 2 channel again shows a positive splitting of the 
order of ~ 50(4) MeV. This suggests that instead of a 
bound state, we appear to be seeing a scattering state in 
the 1 = 2 channel. 

The behaviour observed for the mass differences be- 
tween the 1 = 1 and the two-particle states is illustrated 
in Fig. [7] A positive mass splitting is observed for all 
the five smallest quark masses and increases as the phys- 



FIG. 7: As in Fig. H but for the I(J P ) = 1(1+) state. 



ical regime is approached. The signature of repulsion at 
quark masses near the physical regime would imply no 
evidence of the resonance in the 1 = 1 channel. This sug- 
gests that the 1=1 tetraquark system (in the quenched 
approximation) is more complicated and different from 
the two-hadron system otherwise. Adopting a linear 
chiral extrapolation, the mass of the I(J P ) = 1(1 + ) 
tetraquark state is estimated as M( q =1 = 1360(30) MeV 
which is lower than the mass of the experimentally ob- 
served a (1474). 



TABLE II: Ratios between the mass difference, AM^ q and 



m p at various pion masses. 



AU(GeV) 



AM£°/M P AMj- 2 /M p 



0.822 
0.751 
0.672 
0.616 
0.550 
0.479 



-0.0096(43) 
-0.0016(54) 
-0.0203(71) 
-0.0206(88) 
-0.022(10) 
-0.021(13) 



0.0660(61) 
0.0674(67) 
0.0523(86) 
0.0568 (82) 
0.0650(81) 
0.0664(92) 



0.0054(36) 
0.0278(53) 
0.0436(66) 
0.0562(78) 
0.0669(91) 
0.0731(94) 



As mentioned above, to obtain a definitive result for 
the signature of a possible scalar tetraquark bound state 
with 1 = will require the implementation of the heavy 
quark limit. The heavy quark mass suppresses relativis- 
tic effects, which complicates the interpretation of light- 
quark states. By giving the quarks a larger mass would 
alter threshold, which in turn would affect the manifes- 
tation of the bound state. To confirm this, we calcu- 
late the mass of the lowest Aq state with varying quark 
mass. We allow the quark mass to be larger (hundreds 
of MeV) so that the continuum threshold for the decay 
q 2 q 2 — > {qq){qq) is elevated. The resulting extracted 
masses and mass differences are tabulated in Table IIIII 
and shown in Fig. [8j respectively. 

The ground-state masses in the scalar channel are 
slightly larger than those obtained at the small quark 
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FIG. 8: Chiral extrapolation of the mass ratios at heavy 
quark masses. The solid triangles denote the 1 = 0, and solid 
diamonds represent the 7=1 tetraquark states, respectively. 
The dashed curves are linear fits to the data. 



masses. The mass difference is consistent with or slightly 
below the mass of the lowest two-pion state at the four 
smallest quark masses. The mass difference appears to be 
approaching a positive constant near the physical limit 
and the tetraquark masses are ~ 50 MeV larger than the 
s-wave two-pion state. This suggests that instead of a 
bound state, we appear to be seeing a scattering state in 
the 1 = channel. The lattice data for the 1 = 1 channel 
appear to be well above the increased two-particle state 
threshold and remains approximately constant with mj. 
In the chiral limit, the difference between the I = 1 Aq 
state and the two-pion threshold is ~ 80 MeV. Again, the 
positive mass difference could be a signature of repulsion 
in this channel. We conclude that there are no surprises 
in the 1 = 1 channel -no evidence of attraction and hence 
no indication of a resonance which could be interpreted 
as the 4g state. 



TABLE III: Ratios between the mass difference, AM4, and 
the m p at large quark masses. 



M^GeV) 


AMi=°/M p 


AMi=7M p 


0.947 


0.0006(30) 


0.0429(42) 


0.896 


0.0002(41) 


0.0474(51) 


0.830 


-0.0188(55) 


0.0503(67) 


0.761 


-0.0137(64) 


0.0538(74) 


0.680 


-0.0093(78) 


0.0570(70) 


0.641 


-0.0085(84) 


0.0594(82) 



Finally, we compute the potential Vi q as a function of 
the distance between the two diquark, R, and the inter- 
nal diquark separation D. We calculate the ground-state 
potential from the behaviour of the Wilson loop, (Wi q ), 
at large t for symmetric planar loops with di = d 2 = 
c?3 = g?4 = D. The plateau for representative values of 
R and D arc shown in Fig. [5] As a result of smearing 



the plateaus are seen at earlier t, which indicates a good 
overlap with the state in question. The errors on all our 
data points are jackknife errors. In our analysis, we cal- 
culate the effective 4q potential by choosing a fit which 
has x 2 /d.o.f < 1 and is insensitive to the fit parameters. 



V ef,<t) 



FIG. 9: Effective potential as a function of t in lattice units. 
The graph corresponds to the internal diquark separation D = 
2 and diquark separation R = 4. 



Fig. (|10p shows the tetraquark potential as a function 
of the diquark separation R for different values of the 
internal diquark distance D. At small values of diquark 
separation and R < D, the data seem to coincide and the 
tetraquark potential converges to the two-meson Ansatz. 
Thus for these geometries the system corresponds to the 
two-meson state. For R > D the tetraquark poten- 
tial is lower than 2V qq (R) but seems to be in excellent 
agreement with one-gluon exchange Coulomb plus multi- 
Y Ansatz [23j . Comparing the V^ q and 2V qq , we see that 
the tetraquark potential starts as a sum of the two-meson 
potential and then cross over to approach the connected 
4g state. This would indicate that ground-state config- 
uration is largely a tetraquark state for large R. The 
flip-flop between the Aq state into the two-meson state 
seems nontrivial. 

Now all prerequisites are available to investigate the 
stability of the 4q state. We notice that when qq and 
qq are well separated, all the four quarks are linked by 
connected V-shaped flux tube, where the total flux-tube 
length is minimised. This would mean that the examined 
Aq system can be qualitatively explained as a connected 
4g state. On the other hand, when the nearest quark and 
antiquark are spatially close, the system is described as 
a two-pion state rather than a Aq system. Therefore, this 
type of flip-flop between the Aq and the two-pion state 
around the cross-over can be interpreted as a flux-tube 
recombination between the two mesons. 
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FIG. 10: Tetraquark potential Vi q for symmetric planar con- 
figurations at D = 2 (open circles) and D = 4 (open trian- 
gles), respectively. The solid curve represents the one-gluon 
exchange plus multi-Y Ansatz, and the dashed curve the two- 
meson Ansatz. 

V. SUMMARY AND CONCLUSION 

We have presented the results of our investigation of 
the tetraquark systems in improved anisotropic lattice 
QCD in the quenched approximation on relatively light- 
and heavy-quark masses. Our analysis takes into account 
all possible uncertainties, such as statistical uncertain- 
ties, finite size, and quenching errors when performing 
the chiral extrapolations. The masses of the J = 
and J = 1 states were computed using field operators, 
which are motivated by the 7T7t and diquark structure. 
By analysing the correlation matrix we determined the 
masses of isospin 0, 1 and 2 channels. In the region of 
pion mass which we are able to access, we saw the evi- 
dence of attraction for the 1 = channel. However, at 
relatively large quark masses, the manifestation of the 
bound state changes and 1 = state appears to be the 
scattering state. 



The observation of repulsion at both light- and heavy- 
quark masses is particularly evident in the 1 = 1 and 
1 = 2 states. Both showed a positive splitting between 
tetraquark and two-pion states, suggesting repulsion as 
opposed to attraction. The evidence of repulsion could be 
associated with the absence of a resonance in both these 
channels. Moreover, in the case of I(J P ) = 1(1 + ) and 
0(2 + ), which are more exotic than the 0(0 + ) state, the 
interpolators had sufficient overlap to allow for a success- 
ful correlation matrix analysis and produced no evidence 
of attraction that one might expect if a tetraquark state 
existed. This produces further evidence that there are no 
spatially-localised lowest- lying tetraquark states. 

From our static tetraquark potential analysis we con- 
clude that we showed that at small internal separa- 
tion between diquarks and with internal diquark dis- 
tance larger than the diquark separation (R < D), the 
tetraquark system is a two-meson state and its static po- 
tential is approximately given by the 2V q q potential. For 
R > D, the static potential is described by confining Vi q . 
In summary, these observations suggest that, the 4g sys- 
tem behaves as a multiquark state at diquark separation 
larger than the internal diquark distance. This behaviour 
is expected to hold for general diquark configurations. It 
will be essential to explore unquenched QCD to establish 
the existence or the absence of a tetraquark resonance in 
lattice QCD before one rules out the possible existence 
of a tetraquark state in full QCD. 
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